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Abstract 

Recent high resolution observations by the Chandra X-ray satellite revealed vari- 
ous substructures in hot X-ray emitting plasmas in cluster of galaxies. For example, 
Chandra revealed the existence of sharp discontinuities in the surface brightness at 
the leading edge of subclusters in merging clusters (e.g., Abell 3667), where the tem- 
perature drops sharply across the fronts. These sharp edges are called cold fronts. 
We present results of three-dimensional (3D) magnetohydrodynamic simulations of 
the interaction between a dense subcluster plasma and ambient magnetized intra- 
cluster medium. Anisotropic heat conduction along magnetic field lines is included. 
At the initial state, magnetic fields are assumed to be uniform and transverse to the 
motion of the dense subcluster. Since magnetic fields ahead of the subcluster slip 
toward the third direction in the 3D case, the strength of magnetic fields in this re- 
gion can be reduced compared to that in the 2D case. Nevertheless, a cold front can 
be maintained because the magnetic field lines wrapping around the forehead of the 
subcluster suppress the heat conduction across them. On the other hand, when the 
magnetic field is absent, a cold front cannot be maintained because isotropic heat 
conduction from the hot ambient plasma rapidly heats the cold subcluster plasma. 
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1 Introduction 



In recent years, high spatial resolution observations of clusters of galaxies 
by Chandra revealed t he existence of sharp discontinuities of X-ray in tensity 
called cold fronts (e.g., Markevitch et al. . 20001 : Vikhlinin et al. . 2001b| ). Their 
existence is challenging because in high temperature plasma, high thermal con- 
ductiv ity quickly smears such a discontinuity in temperature. lEttori and Fabian 
(2000) pointed out that heat conductivity across the front should be reduced 



from the classical Spitzer value, ks p = 5 x 10 7 T 5 / 2 ergs 1 cm 1 K 1 (jSpitzer, 
by several orders of magnitude. Magnetic fields enable such suppression. 
Cold fronts ma y also b e subject to the Kelvin- Helmholtz (K-H) instability. 



Vikhlinin et al.l ([2001aU bl) suggested that turbulent magnetic fields stretched 
along the cold front will reduce the growth rate of the K-H instability at the 
interface between a subcluster plasma and an ambient plasma. 

The magnetic fie ld strength in the core of clusters of ga l axies can be the 
order of /iG (e.g., Kronbere . 1994 : Carilli and Tavlor . 20021 ). Johnston-Hollitt 



(12004 ) summarized the observations of the field strength in several regions of 
Abell 3667 by using different techniques. They showed that its strength in the 
cluster core is about 1 — 2 /iG, and that magnetic fields are tangled on scales 
of roughly 100 kpc. 



Let us estimate the heat conduction time scale in cluster of galaxies. It has 
been sugge sted that the heat conductio n in cluster of galaxies can be very effi- 
cient (e.g.. lTakahara arid Ikeuchil . ll977l ) because the Coulomb mean free path, 



l c ~ 5 (kT/5 keV) 2 (n/10 3 cm 3 ) 1 kpc, is large. Here, n is the number density. 
The time required for heat to diffuse by conduction across a length L is given 
by T Sp ~ pL 2 /n Sp ~ 3 x 10 7 (kT/5 keV)- 5 / 2 (n/10" 3 cm- 3 )(L/100 kpc) 2 yr, 
where p is the density and Ks P is the Spitzer conductivity. Since the thick- 
ness of the cold front is much smaller than 100 kpc, they will disappear in a 
time scale shorter than 10 7 yr. If magnetic fields exist, however, the charac- 
teristic length of the heat exchange across magnetic field lines is reduced to 
the Larmor radius, = 2500 (B/l fiG)~~ 1 (kT/5 keV) 1 / 2 km. Thus, the heat 
conduction across the magnetic field is almost suppressed. 



Earlier numerical studies of cold fronts ( e.g., Heinz et al. , 2003| ) did not in- 
clude magnetic fields and heat conduction. lAsai et aL ( 2004f ) first reported the 
results of two-dimensional (2D) magnetohydrodynamic (MHD) simulations of 
cold fronts including both magnetic fields and anisotropic heat conduction. 
They showed that magnetic fields are stretched and elongated ahead of the 
front. Since these magnetic fields suppress the heat conduction across them, 
the contact discontinuity (a cold front) between the cold subcluster plasma and 
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the hot ambient plasma can be maintained. The cold subcluster is wrapped 
by magnetic fields and protected from being heated by the heat conduction. 

In this paper, we focus on the three-dimensional (3D) effects. In the 3D case, 
since magnetic fields compressed in front of the cold front can expand in the 
third direction or slip out the front surface, the strength of the magnetic fields 
in this region may be reduced compared to that in the 2D case. We would 
like to show that even when we include these 3D effects, cold fronts can be 
maintained. 



2 Simulation model 



We simulated the time evolution of a cluster plasma in a frame comoving with 
the subcluster. The basic equations are as follows: 



dp 
di 

P 

OB 

~dt 



+ V -(pv) = 0, 
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where p, v, p, B, and ip are the density, velocity, pressure, magnetic fields, and 
gravitational potential, respectively. We use the specific heat ratio 7 = 5/3. 
The subscript || denotes the components parallel to the magnetic field lines. 
We assume that heat is conducted only along the field lines. We solved ideal 
MHD equations in a Cartesian coordinate system (2, y, z) by a modified Lax- 
Wendroff method with artificial viscosity. The heat conduction term in the 
energy equation is s olved by the implicit red and black successive overre- 



laxation method (see lYokovama and Shibatal . l200ll for detail). The radiative 



cooling term is not included. The units of length, velocity, density, pressure, 
temperature, and time in our simulations are r$ = 250 kpc, vq = 800 km s -1 , 
p = 5 x 10~ 27 gcm -3 , p = 3 x 10~ n erg cm -3 , kT = 4keV, and t = r /v = 
3 x 10 8 yr, respectively. 
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Figure 1 shows the initial density distribution. Solid lines and arrows show 
magnetic field lines and velocity vectors in the y = plane, respectively. We 
assume that a spherical isothermal low-temperature (kT in = 4keV) plasma is 
confined by the gravitational potential of a subcluster. The low-temperature 
plasma is embedded in the less-dense (p out = 0.25 po) uniform, hot (kT out = 
8keV) plasma. Here the subscripts "in" and "out" denote the values inside 
and outside the subcluster, respectively. 



We assume tha t the density distribution of the su bcluster is given by the 
/5-model profile ( Cavaliere and Fusco-FemianoL 197fi| ). 







1 + 









-3/3/2 



(5) 



where r = (x 2 + y 2 + z 2 ) 1 ^ 2 , (3 = 2/3, the maximum density p c = 10~ 26 gcm -3 , 
and the core radius is r c = 290 kpc. The subcluster is initially in hydrostatic 
equilibrium under the gravitational potential fixed throughout the simulation. 

We assume that ambient plasma initially has a uniform speed with Mach 
number M = v x /c sont = 1, where c sout is the ambient sound speed. The 
Mach number with respect to the sound velocity inside the subcluster is M' = 
v x /c sin = y/2. 

Table 1 shows model parameters. An important parameter is the plasma beta 
(/3 ) defined as the ratio of the ambient gas pressure to the magnetic pressure. 
When p = Pgas/Pmag = 100, the magnetic field strength is B ~ 1.5 /zG. 
Models HC2 (2D) and HC3 (3D) are non-magnetic models with isotropic heat 
conduction, and model MC2a (2D) and MC3a (3D) are models with magnetic 
fields (/?o = 100) and anisotropic heat conduction. Model MC2b (2D) is a 
model with weak magnetic fields (/?o = 1000, B ~ 0.5 /xG) and anisotropic 
heat conduction. Model MC3b (3D) is a model with weaker magnetic fields 
(/3o = 10 3 , B ~ 0.15 fiG). We note that the initial magnetic fields are parallel 
to the 2-axis and perpendicular to the motion of the subcluster. 

When magnetic fields exist, heat is conducted only in the direction parallel to 
the field lines. Along the magnetic fields, the heat conductivity is taken to be 
that of the Spitzer conductivity, K |l = 5 x 10~ 7 T 5 / 2 erg s 1 cm" ^K -1 . Mean- 
while the conductivity across the field lines is n± = 0. The conduction time 
scale along the field lines is r K ~ pL 2 / 'k\\ ~ 3x 10 7 (kT/5 keV)~ 5 / 2 (n/10~ 3 cm -3 ) 
(L/100kpc) 2 yr. 

For boundary conditions, the left boundary at x = —5 is taken to be a fixed 
boundary, and other boundaries are taken to be free boundaries where waves 
can be transmitted. 
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Fig. 1. Initial density distribution visualized as logarithm of the density. 
Table 1 

Models and parameters. Heat conduction is isotropic in models HC2 and HC3. In 
other models, heat conducts only along magnetic fields. 



Model 


A) 


K 


A;T in [keV] 


A;T out [keV] simulation box 


number of grid points 


HC2 


oo 


K 


4 


8 (2.5Mpc) 2 


1200 2 


MC2a 


100 


K|| 


4 


8 (2.5Mpc) 2 


1200 2 


MC2b 


1000 


K|| 


4 


8 (2.5Mpc) 2 


1200 2 


HC3 


oo 


K 


4 


8 (2.5Mpc) 3 


280 3 


MC3a 


100 


«|| 


4 


8 (2.5Mpc) 3 


280 3 


MC3b 


10 4 


K„ 


4 


8 (2.5Mpc) 3 


280 3 


3 Results 



3. 1 Configuration of magnetic fields 



Figure 2 shows the magnetic fields for model MC3a at t = 1 Gyr. Solid curves 
and color contours show the magnetic field lines and density distribution, 
respectively. The subcluster plasma is wrapped by the magnetic field lines. 
Ahead of the subcluster, the magnetic field lines are almost parallel to the 
contact surface between the subcluster plasma and the ambient plasma. We 
can identify magnetic field lines slipping in the y-direction. 

Figure 3 shows isocontours of magnetic field strength in the y = plane for 
models MC3a (left, 3D) and MC2a (right, 2D) at t = 1 Gyr, respectively. 
Arrows show velocity vectors. We note that the initial magnetic field lines 
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Fig. 2. Magnetic field lines at t = 1 Gyr. Color contours show the logarithmic density 
distribution. The solid curves show magnetic field lines. 

are parallel to the z-axis. The magnetic field strength is enhanced ahead of 
the subcluster via compression. Compared to the 2D model (right), magnetic 
fields in the 3D model are weaker ahead of the subcluster but their strength 
is enhanced at the tail of the subcluster because the plasma flows advect the 
magnetic fields toward the tail and they pile up in the tail. 



The left panel of Figure 4 shows the distribution of magnetic field strength 
at t = 1 Gyr along the x-axis. The solid curve shows their strength for model 
MC3a (3D, (3q = 100). The dotted line and the long-dashed line show their 
strength for model MC2a (2D, f3 = 100) and MC2b (2D, f3 = 1000), respec- 
tively. Their strength in front of the cold front is enhanced more than a factor 
3 from the initial state. The right panel in Figure 4 shows the time evolution 
of the total magnetic energy integrated in the whole simulation region. The 
magnetic energy increases gradually due to magnetic field compression around 
the subcluster. 



3.2 Effects of magnetic fields on heat conduction 



Let us compare the results for model MC3a (with heat conduction) and model 
HC3 (without heat conduction). We show in Figure 5 distributions of the pro- 
jected density, temperature, and X-ray intensity obtained by integrating these 
quantities in the ^/-direction. X-ray intensity are visualized from simulation 
results as logarithm of the thermal bremsstrahlung emissivity, ~ p 2 . In model 
HC3 (lower panels), heat is conducted isotropi cally. Since the conducti on time 
scale across the cold front with width ~ 5kpc ( Vikhlinin et ah . 2001b| ) is very 
short (r K ~ 10 7 yr), the cool component inside the subcluster is heated by the 
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Magnetic field strength 

/ = 3.0/(j B f = 3.O( B 




model MC3a x model MC2a * 



Fig. 3. Comparison of the magnetic field strength at t = 1 Gyr between 3D (left) 
and 2D (right). The left and right panels show the result for models MC3a and 
MC2a at y = plane, respectively. Arrows show the velocity vectors. 




Fig. 4. Left: Distributions of magnetic field strength at t = 1 Gyr. The solid curve 
shows magnetic field strength for model MC3a (3D, (5q = 100) along the x-axis. The 
dotted curve and the long-dashed curve show the strength for model MC2a (2D, 
/9b = 100) and MC2b (2D, (3 = 1000) along the x-axis, respectively. The dashed 
line and dash-dotted line show the initial magnetic field strength for models with 
/3q = 100 and (3q = 1000, respectively. Right: Time evolution of the magnetic energy 
(Emag) integrated in the whole simulation region. The increase with respect to the 
initial kinetic energy, A£ mag /£ kin (0) = (E ma , g (t) - E mag (0)) / 'E kin (0) , is plotted. The 
solid line, dotted line, and long-dashed line are for models MC3a (3D, (3q = 100), 
MC2a (2D, /3 = 100), and MC2b (2D, (3 = 1000), respectively. 

hot ambient plasma and evaporates rapidly. Meanwhile in model MC3a (top 
panels), heat is conducted only in the direction parallel to the field lines. As 
we showed in §3.1, since the subcluster is wrapped by magnetic field lines, the 
cool component inside the subcluster is protected from heating. We can see 
this protection remarkably in the temperature distributions (middle panels). 
Thus the discontinuity remains sharp. 



7 



DltlsiI.v 



I cjii ipcrriLurc 



-=1 -2 
model HC3 I 



1 "1 


1 * 




7.1 


LOO. 


K (1 


50.0 


-2- 
-1- 



-1 -2 I) ?. i 

model ML'jn 
j 3.0 <a 



150.1 



50.0 I 



-4 -2 2 
model _\lC3a - 1 

f 3.0 ii) 



X-ray intensity 
kT r=3.()f t | 




-A -2 2 
uiodcL ML'jn 1 




• 



200. n 



X-FH}' 

.2(1(1. 



-■1 J (I ?. -1 
inodcl HC3 * 



-l -7. I) l 
model HCj « 



Fig. 5. Distribution of projected density (left), temperature (middle), and X-ray 
intensity (right) integrated at t = lGyr along the y-axis. The top panels show results 
for model MC3a (with magnetic fields and heat conduction). The bottom panels 
show results for model HC3 (without magnetic fields and with heat conduction). 



Figure 6 shows the temperature profiles at t = lGyr along the x-axis. The 
left panel shows those for models MC3a (solid line) and HC3 (dotted line). 
The right panel shows those for models MC2a (solid line), MC2b (dash-dotted 
line), and HC2 (dotted line). In both panels, the dashed line shows the initial 
temperature profile. It is obvious that magnetic fields along the front enables 
the cold front to exist for lGyr. 



We can also identify shock waves. Numerical results reproduced the bow shock 
produced by supersonic motion upstream of the cold front observed in A3667 
(see Figure 5). The interaction of the subcluster plasma with the ambient 
plasma generates another shock propageting inside the subcluster. This shock 
exists at x = 2.5 in models MC3a and MC2a (solid lines) in Figure 6, and 
also visible in the left panel of Figure 4. These shocks drive the adiabatic 
expansion of the subcluster and cools t he subcluster plas ma. Such behaviors 
are consistent with those in 2D models ijAsai et all 120041 ) . 
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Fig. 6. Comparison of the temperature distributions between 3D (left) and 2D 
(right) at t = 1 Gyr along the x-axis. The left panel shows the results for mod- 
els MC3a (solid line), MC3b (dash-dotted line), and HC3 (dotted line). The right 
panel shows the results for models MC2a (solid line), MC2b (dash-dotted line), and 
HC2 (dotted line). The dashed lines in both panels show the initial temperature 
distribution. 

4 Discussion 



4-1 Importance of magnetic fields for the existence of cold fronts 



We have shown that magnetic fields enable the existence of cold fronts by sup- 
pressing the heat conduction a cross the front. Previous hydrodynamic simu- 
lations (e.g., Heinz et al. . 2003) could reproduce the cold fronts because they 
neglected the heat conduction. When heat conduction is included, cold fronts 
disappear in time scales of ~ 10 8 yr unless the heat conduction is suppressed 
as we showed in models HC2 and HC3 (see Figures 5 and 6). Magnetic fields 
much smaller than that in model MC3a can suppress heat conduction. We 
confirmed that weak magnetic fields for model MC3b (flo — 10 4 , B ~ 0.15 /xG) 
can suppress the heat conduction across the front as long as they cover the 
front. Even when large-scale coherent magnetic fields do not exist, turbulent 
magnetic fields can suppress the heat conduction ac ross the front; they are 
stretched and elo ngated parallel to the front (see also IVikhlinin et all 12001a ; 



Asai et all 12004 b 



Vikhlinin et al. 



(|2001al lbh suggested that magnetic fields are important for 
the existen ce of cold f r onts b ecause they can suppress the growth of the K-H 
instability. I Asai et ajj (J2004J) showed that in models without heat conduction 
and without magnetic fields, the growth of the K-H instability is not prominent 
and does not disrupt the cold front ahead of the moving subcluster. In the 
sides of the subcluster, the K-H instability creates eddies which detach and 
propagate downstream. Magnetic fields can suppress the formation of such 
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eddies. 



4-2 Possibility of subcluster disruption by magnetic fields accumulating ahead 
of the subcluster 



Let us now disc uss the possibil i ty of t he disruption of the subcluster plasma by 
magnetic fields. Gregori et al. ( 2000h presented the results of 3D simulations 
of moderate supersonic cloud motion in a magnetized interstellar medium. 
They assumed that the cloud has uniform density. They showed that magnetic 
fields accumulating ahead of the cloud triggers the Rayleigh- Taylor (R-T) 
instability and disrupts the cloud. Our numerical results, however, does not 
show the growth of such R-T instability. These different behaviors are due 
to a difference of the density distribution of the cloud. Since the subcluster 
plasma in our simulations has a dense core because it is confined by the gravity, 
magnetic fields ahead of the clump can not deform the density distribution 
inside the subcluster. 



4-3 Three-dimensional effects 



We have shown that in the 3D case, since magnetic fields ahead of the sub- 
cluster slip along the contact surface between the subcluster plasma and the 
ambient plasma, the amplification of their strength is smaller than that in the 
2D case. This weaker magnetic field, however, is sufficient to suppress the heat 
conduction across the front. 

Since the area of the contact surface is larger in the 3D case than that in 
the 2D case. Thus cold subcluster plasma can be heated more efficiently than 
in the 2D case. We confirmed this by comparing results for model HC2 and 
HC3 (Figure 6). Compared with model HC2 (2D, non- magnetic fields), the 
temperature profile for model HC3 (3D, non-magnetic fields) shows that the 
cold subcluster is heated up earlier. When magnetic fields exist, however, the 
difference becomes less prominent. 
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